ASYMPTOTIC EIGENVALUE DISTRIBUTIONS 
OF BLOCK-TRANSPOSED WISHART MATRICES 

TEODOR BANICA AND ION NECHITA 

Abstract. We study the partial transposition = (id(g)t)VF G Md„(C) of a Wishart 
matrix W € Mrf„(C) of parameters {dn,dm). Our main result is that, with d — )■ oo, 
the law of mW^ is a free difference of free Poisson laws of parameters m{n ± l)/2. 
Motivated by questions in quantum information theory, we also derive necessary and 
sufficient conditions for these measures to be supported on the positive half line. 



Introduction 

The partial transposition oi a. d x d block matrix W G M(i(C) ® M„(C) is the matrix 
obtained by transposing each of the n x n blocks of W. That is, we let = 
(id t)W, where id is the identity of Md(C), and t is the transposition of M„(C). The 
partial transposition operation can be defined using coordinates, as follows. A particular 
decomposition C^" = ^ induces a decomposition Mrf„(C) = Mrf(C) ® M„(C). The 
entries of a matrix W G Mrf(C) (g) M„(C) can be indexed by four indices, i,j G {1, . . . , d} 
identifying the block of A the matrix entry belongs to, and a,b G {1, . . . , n} fixing the 
position of the matrix entry inside the block. Then, the partial transposed matrix 
has coordinates 

Motivated by questions in quantum information theory, the partial transposition op- 
eration for Wishart matrices was studied by Aubrun who showed that, for certain 
special values of the parameters, the empirical spectral distribution of converges in 
moments to a non-centered semicircular distribution. In this paper we discuss the general 
case, our main result being as follows. 

Theorem A. Let W be a complex Wishart matrix of parameters {dn,dm). Then, with 
d ^ oo, the empirical spectral distribution of mW^ converges in moments to a free dif- 
ference of free Poisson distributions of respective parameters m{n ± l)/2. 

Observe that if one applies the transposition map on the first factor of the tensor 
product Mrf(C)(8)Mfi(C) instead of the second one, the spectral distribution of remains 
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unchanged. Indeed, this follows from (t ® id) = (t ® t) o (id ® t) and from the fact that 
applying the global transposition leaves the spectral distribution unchanged. 

The above theorem basically generalizes Aubrun's result in [1], modulo some standard 
free probability facts, regarding the relationship between the free Poisson laws and the 
semicircle laws. Also, as a remark of algebraic interest, in the formal limit m = t/n — )■ 0, 
the law that we obtain is a free difference of free Poisson laws having the same parameter, 
namely t/2, so it is a (modified) free Bessel law in the sense of [2]. 

Our second result regards the various properties of the limiting measure. By using tools 
from free probability theory and from classical analysis, we obtain the following result. 

Theorem B. The limiting measure found in Theorem A has the following properties: 

(1) It has at most one atom, at 0, of mass max{l — mn, 0}. 

(2) It has positive support iff n < m/4 + 1/m and m >2. 

The main motivation for the above results comes from quantum information theory. In 
quantum information theory, the partial transposition map is known to be an entangle- 
ment witness: it allows to test if a quantum state (represented by a positive, unit trace ma- 
trix) is entangled, in the following sense. If a bi-partite quantum state p G Mrf(C)®M„(C) 
is separable (i.e. it can be written as a convex combination of product states pf^ ® pP^), 
then its partial transposition p^ is also a quantum state. However, if p is entangled, then 
p^ may fail to be positive. In the case where p^ is a positive matrix, the quantum state p is 
said to be PPT (Positive Partial Transpose). Hence, separable states are always PPT and 
non-PPT states are necessarily entangled. The equivalence of entanglement and non-PPT 
is known to hold only for total dimension smaller than 6 (2 x 2 or 2 x 3 product systems) 
and it fails for larger dimensions, in the sense that there exist PPT entangled states. In 
the same spirit as in pLj, the results regarding the positivity of the support of the limit 
measure can be interpreted as results about typicality of PPT states for large quantum 
systems. Wishart matrices (normalized to have unit trace) are known to be physically 
reasonable models for random quantum states on a tensor product C^^C", the parameter 
m of the Wishart distribution being related to the size of some environment C^"* needed 
to define the state. So, as a conclusion. Theorem A and Theorem B above indicate that, 
when m > 2 and n < m/A + 1/m, a typical state in (g) C" is PPT. 

The paper is organized as follows: in [1] we present a combinatorial formula for the 
asymptotic moments of mW^, in |2] we find an equation for the corresponding moment 
generating function, and in [3] we state and prove the main result. In sections H] and |5] 
we discuss some of the properties of the asymptotic eigenvalue distributions (examples, 
atoms, density, positivity), and the final section [H] contains a few concluding remarks. 

Acknowledgments. We would like to thank the ANR projects Galoisint and Granma 
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work was done when I.N. was a postdoctoral fellow at the University of Ottawa where he 
was supported by NSERC discovery grants and an ERA. 
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1. Formula for moments 



We recall that a complex Wisliart matrix W of parameters {dn, dm) is defined as 
W = {dm)~^GG* , where G is a x dm matrix with i.i.d. complex Gaussian A/'(0, 1) 
entries. Consider the block-transpose matrix , constructed in the introduction. We 
denote by tr : Mrf„(C) — )■ C the normalized trace, tr(iy) = {dn)~^ ^ Wu. 

We recall a number of results from the combinatorial theory of noncrossing partitions; 
see [7] for a detailed presentation of the theory. For a permutation a G Sp, we introduce 
the following standard notation: 

• #0" is the number of cycles of a; 

• \a\ is its length, defined as the minimal number k such that a can be written as 
a product of k transpositions. The function {(t,tt) — t- |o"~^7r| defines a distance on 
Sp. One has #cr + |cr| = p. 

• e((j) is the number of blocks of even size of a. 

Let 7 G S'p be a fixed arbitrary full cycle (an element of order p in Sp). The set of 
permutations a G Sp which saturate the triangular inequality |cr| + |(T^^7| = I7I = p — 1 
is in bijection with the set NC{p) of noncrossing partitions of {1, . . . ,p}. We call such 
permutations geodesic and we shall not distinguish between a non crossing partition and 
its associated geodesic permutation. We also recall a well known bijection between NC{p) 
and the set NC2{2p) of noncrossing pairings of 2p elements. To a noncrossing partition tt G 
NC{p) we associate an element fat(7r) G NC2{2p) as follows: for each block {ii, ^2, • • • , ik} 
of TT, we add the pairings {2ii — 1, 2ik}, {2ii, 2x2 — 1}, {2^2, 2^3 — 1}, ... , {2zfc_i, 2ik — 1} to 
fat(7r). The inverse operation is given by collapsing the elements 2i — 1, 2z G {!,..., 2p} 
to a single element i G {1, . . . The number of blocks of vr and fat(7r) are related 
by #7r = 7^(fat(7r) V P12), where V is the join operation on the lattice NC2{2p) and 
P12 = (12)(34) ■ ■ ■ (2p — I2p) is the fattened identity permutation. Similarly, #(717) = 
^(fat(7r) V P14), where pi4 is the pairing that corresponds to the fattening of the inverse 
full cycle 7^^(i) = i — 1. More precisely, pi4 pairs an element 2i with 2{i — 1) — 1 = 2i — 3, 
or, equivalently, an element i G {1, . . . , 2p} with i + (— 1)*"'"^3. 

The following combinatorial lemma is essential in the proof of the moment formula in 
Theorem II. 2[ 

Lemma 1.1. For every noncrossing partition tt G NC{p), 



Proof. We use a recurrence over the number of blocks of vr. If vr has just one block, its 
associated geodesic permutation is 7 and one has 



l + e(7r) = #(7r7). 




if p is odd; 
if p is even. 
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Figure 1. The partition fat(7r) V pi4 (top line) has one more block than 
fat(7r') V p']^4 (bottom) when k is even. 



For partitions tt with more than one block, we can assume without loss of generality that 
vr = Ifc U it', where 1^ is a contiguous block of size k. Recall that the number of blocks of 
the (not necessarily geodesic) permutation is given by #(vr7) = #(fat(7r) Vpi4), where 
Pi4 e 112 (2p) is the pair partition which pairs an element i with i + (— 1)*+^3. 

If k is an even number k = 2r, then the partition fat(l2r U tt') V pi4 contains the block 
(14 5 8 ■ ■ ■ 4r — 3 4r) along with the blocks coming from elements of the form 4? + 2, 4i + 3 
from {l,...,4r} and from tt'. Since we are interested only in the number of blocks 
statistic (and not in the contents of each block), we can count the blocks of the join 
of two partitions by drawing them one beneath the other and counting the number of 
connected components of the curve, without taking into account the possible crossings, 
see Figure [H We conclude that #(fat(7r) V P14) = 1 + #(fat(7r') V p'l^) where p'14 is pi4 
restricted to the set {2k + 1, 2fc + 2 . . . , 2p}. If k is odd, = 2r + 1, there is no extra 
block appearing, so #(fat(7r) V P14) = #(fat(7r') V p[^), see Figure [2J □ 

Theorem 1.2. For any p > 1 we have 

lim (E o tr) {mW^Y = V m#^n"(^) 

TreNC{p) 

where #(•) and e(-) denote the number of blocks and the number of blocks of even size 
statistics. 
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Figure 2. The partition fat(7r) V pu (left) has the same number of blocks 
as fat(7r') V (right) when k is odd. 

Proof. The matrix elements of the partial transpose matrix are given by: 

d m 
k=l c=l 

This gives: 



isClsyts+lO-a+l 

ii,...,ip=l ai,...,ap=l s=l 

d n p 



ii,...,ip=l ai,...,ap=l s=l 

d n p d m 

^{dn)-\dm)-^ E E n E E G,^a.^ujsbA.^.a.,.,.. 
ii,...,ip=l ai,...,ap=l s=l ji,...,jp=l bi,...,hp=\ 

After interchanging the product with the last two sums, the average of the general term 
can be computed by the Wick rule, namely: 

(p \ p 

Let 7 G 5*^ be the full cycle 7=(12...p)~^. The general factor in the above product 
is 1 if and only if the following four conditions are simultaneously satisfied: 

7~''^7r < keri, 7r7<kera, tt < ker j, 7r<ker6, 

where the kernel of a function / : {1, . . . ,p} — )■ M is the partition ker/ where i and j are 
in the same block iff f{i) = f{j). Counting the number of free parameters in the above 
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equation, we obtain: 

The exponent of d in the last expression on the right is 

iV(7r) = #7r + #(7-V)-p-l 
= p - 1 - (|7r| + |7~V|) 
= p - 1 - (|7r| + |7r"-^7|) 

and, as explained in the beginning of this section, this quantity is known to be < 0, with 
equality iff n is geodesic, hence associated to a noncrossing partition. We get: 

(E o tT)[{W^Y] = (1 + 0{d-^))m-Pn-^ ^ m#"n#("^). 

7reiVC(p) 

By Lemma [1.1^ we have #(vr7) = e(7r) + 1, which gives the result. □ 

2. Moment generating fungtion 

In this section we find an equation for the moment generating function of the asymptotic 
law of mW^ . This moment generating function is defined by: 

F(z) = lim (E o tr) ( ^ -] . 

Equivalently, F{z) = Yl'^=o ^p^^i where Mp is the moment computed in Theorem 11.21 
Note that one can crudely upper bound the moment Mp by A^{mnY so that the moment 
generating function above is defined at least for z small engouh, \z\ < (4mn)~^. 

Theorem 2.1. The moment generating function of mW^ satisfies the equation: 

(F - 1)(1 - z^F'^) = mzF{l + nzF) 

Proof. We use Theorem 11.21 If we denote by N{p, b, e) the number of partitions in NC{p) 
having b blocks and e even blocks, we have: 

oo oo oo oo 

F = l + J2 Yl zPm*^n''^^^ = 1 + 1] 1] 5Z z''m^n^N{p, b, e) . 

p=l it€NC{p) P=1 b=0 e=0 
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Let us try to find a recurrence formula for the numbers N{p,b,e). If we look at the 
block containing 1, this block must have r > other legs, and we get: 

N{p,b,e) = Y, Yl Y 5^ Ar(pi,6i,ei)...iV(p,+ 

re2N p=Spi+r+l b=T,bi+l e=Sei 

+ Y Y Y Y N{pi,bi,ei) . . .N{p 

re2N+l p=Spi+r+l b=T.b^+l e=Ee,+l 

Here Pi, ■ ■ ■ ,Pr+i are the number of points between the legs of the block containing 1, 
so that we have p = (pi + . . . + Pr+i) + r + 1, and the whole sum is split over two cases (r 
even or odd), because the parity of r affects the number of even blocks of our partition. 

Now by multiplying everything by a z^m^n'^ factor, and by carefully distributing the 
various powers of z, m, b on the right, we obtain the following formula: 

r+1 

zPm^7fN{p, b,e)=mY ^'^^ Y Y YX\ z^^m^^7f^N{pi, k, a) 

re2N p=Spi+r+l fe=Sbj+l e=Eei i=l 

r+1 

+ mn Y ^'^^ Y Y Y X\zP^'m^'n^'N{piM,ei). 

re2N+l p=Spi+r+l b=Efei+l e=Sej+l i=l 

Let us sum now all these equalities, over all p > 1 and over all 6, e > 0. According to 
the definition of F, at left we obtain F — 1. As for the two sums appearing on the right 
(i.e., at right of the two z"^^^ factors), when summing them over all p > 1 and over all 
6, e > 0, we obtain in both cases F"^^^. So, we have the following formula: 

F-l=mY {zPy^^ + mn Y i^FY^^ 

re2N rG2N+l 

zF z^F^ 

= m — - + mn — - 

1 - z^F^ 1 - z^F^ 

^ 1 + nzF 

This gives the formula in the statement. □ 

3. Limiting distribution and free Poisson laws 

In this section we state and prove our main result. Let us first recall that for a complex 
Wishart matrix W of parameters {dn,dm), the eigenvalue distribution of tW, with t = 
m/n, converges in moments with c? — oo to the Marchenko-Pastur law [6], given by 



J4t-{x-i- ty 

TTf = max(l — t, 0)Oq + — 1 

where a = {y/i - 1)^ and b = {Vt + if. 



-ia,fe](a;) (ia;. 
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In order to deal with the block-transposed Wishart matrices W^, we will need a free 
probability point of view on tt^. So, let us recall from [8] that a noncommutative proba- 
bility space is a pair {A, (p), where yl is a C*-algebra, and : A — )■ C is a positive unital 
linear form. The law of a self-adjoint element a & A is the probability measure on the 
spectrum of a (which is a compact subset of M) having as moments the numbers ip{a^). 

Two subalgebras B,C & A are called free if ip{. . . hiCihi^iCi^i . . .) = whenever h E B 
and Q G C satisfy '^{hi) = if{ci) = 0. Two elements b,c E A are called free whenever 
the algebras B =< b > and C =< c > that they generate are free. Finally, the free 
convolution operation ffl is defined as follows: if /i, u are compactly supported probabihty 
measures on M, then /x ffl z/ is the law of 6 + c, where b, c are free, having laws /i, u, see [S]. 

One of the remarkable results of the theory is that when performing a "free Poisson 
limit" , we obtain precisely the Marchenko-Pastur law. That is, we have: 



For this reason, the Marchenko-Pastur law is also called "free Poisson law" [7]. 

In what follows we will find a free probability interpretation of our main result so far, 
namely Theorem 12.11 We will prove that mW^ has the same law as a — 6, where a, b are 
free variables following the laws vr^, vr^, with s = m{n + l)/2 and t = m{n — l)/2. 

Theorem 3.1. With d — )■ oo, the empirical spectral distribution of mW^ converges in 
moments to a probability measure fim,n, which is the free difference of free Poisson distri- 
butions of parameters m{n ± l)/2. 

Proof. We use Voiculescu's _R-transform, which is defined as follows [8j. Let be a 
compactly supported probability measure on M, having moment generating function F. 
The Cauchy transform of /i is then G(^) = i^^F{^~^). We denote by K the formal inverse 
of G, given by G{K{z)) = K{G{z)) = z, and then we write K{z) = R{z) + z~^. The 
function R = R^, called /^-transform of fi, has the remarkable property R^wu = R^i + Ru- 
in other words, R is the free analogue of the logarithm of the Fourier transform. 
Consider now the equation of F, found in Theorem 12.11 




(F - 1)(1 - z^F^) = mzF{l + nzF). 
With z and F (so that zF G) we obtain: 

i^G - 1)(1 - G^) = mG{l + nG). 
Now with C, K and G — ?■ z we obtain: 

{zK - 1)(1 - z'^) = mz(l + nz). 
Finally, with K ^ R + z~^ we obtain: 

zR{l — z^) = mz{l + nz). 
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Thus the i?-transform of the asymptotic law of mW^ 

1 + nz m f n + 1 i 




is given by: 
n-l\ 



l + zj' 



Let now a, b be free variables, following free Poisson laws of parameters s, t. We have 
Ra = — z) and Rh = t/{l — 2), and by using the general dilation formula Rqb{z) = 
qRh{qz) at g = — 1 (see j7]) we deduce that we have R^h = + z). So, we have: 



Now since with s = m{n + l)/2 and t = m{n — l)/2 we obtain the above formula of R, 
we are done. □ 



In the reminder of this paper we study the limiting measures found in Theorem 13.11 It 
is convenient to enlarge our study to all the possible values of m, n, as follows. 

Definition 4.1. To any real numbers m > and n > 1 we associate the measure fim,n 
which is the free difference of free Poisson distributions of parameters m{n ± l)/2. 

As a first remark, the various formulae found in this paper (e.g. moment formula in 
Theorem II. 2[ and various equations for F, G, K, R in the proof of Theorem 13. ip are of 
valid in this more general setting. Let us collect some useful formulae here. 

Proposition 4.2. The measure fim,n has the following properties: 

(1) It has mean m and variance mn. 

(2) Its Cauchy transform satisfies {C,G — 1)(1 — G^) = mG{l + nG). 

(3) Its R-transform is given by R{z) = m{l + nz)/{l — z"^). 

Proof. The first statement follows from the fact that both the average and the variance 
of a free Poisson law of parameter t are t. The last two items follow from the proof of 
Theorem 13. II above, by proceeding backwards. □ 

Regarding now the main particular cases of measures of type /im,n, corresponding to 
some previously known random matrix computations, the situation is as follows. 

Theorem 4.3. The measures ^m,n o,^^ o,s follows: 

(1) At n = 1 we have the Marchenko-Pastur laws: fim,i = T^m- 

(2) With m = an 00 we obtain Aubrun's shifted semicircles [1] . 

(3) With m = t/n we obtain a modified free Bessel distribution Ti2,t- 

Proof. All the results are clear either from Definition 14. or from Proposition 14.21 
(1) This follows from Definition 14. because at n = 1 we have m{n ± l)/2 = m, 0. 



4. Examples of limiting measures 
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(2) The precise statement here is that for any P = a ^ > we have the formula 



hm Dm-i{flrn,/3m) = SC{1,(3) 
m—^co 

where D is the dilation operation (if fi is the law of a random variable a then Dq{fi) is 
by definition the law of the random variable qa) and where SC{1, /3) is the semicircle law 
centered at 1, having support [1 — 2a/5, 1 + 2a/^. 

Indeed, let F = r(^) be the Cauchy transform of u = Dm-i(/^m,n)- The Cauchy trans- 
form of /im,n = Drn{v) is then given by G{m^) = m'^T, and by making the replacements 
^ — 7- mC, and G — ?■ m^^T in the equation of G in Proposition 14. 2^ we obtain: 

(^r - i)(i - m-^r^) = r(i + nm-^r). 

In the limit n = (3m — )■ oo, as indicated above, this equation becomes: 

- 1 = r(i + (3T) =^ /3r2 + (1 - ^)r + 1 = 

^-l±^(l-02_4^ 

2/3 



By applying now the Stieltjes inversion formula, we get the following density: 

fix) 



- (1 - 



X] 



2(371 



But this is exactly the density of the semicircle law SC{1, (3), and we are done. 

(3) This follows directly from Definition 14.11 because with m = t/n — )■ Owe have 
m{n ± l)/2 = t/2,t/2, and these parameters are the same as those for the modified free 
Bessel law 7r2,j, see p]. Definition 7.2. □ 

Observe that the formula (1) above is in tune with the Marchenko-Pastur theorem, as 
formulated in section |3l indeed, at n = 1 the partial transposition operation is trivial, so 
Theorem 13. II computes the asymptotic law of the rescaled Wishart matrix mW. 

Also, (2) above agrees of course with Aubrun's computation in [1], with Theorem 13.11 
generalizing Aubrun's result. This result can be seen as well to follow from the following 
heuristic argument using the free central limit theorem. Using the free Poisson limit 
theorem, one can write (we denote a Bernoulli random variable by b{p, x) = {l—p)6o+pSx) 



lim h 

k—^oo 



171(71 



6i ffl lim 

A;— >oo 



2k rr 

m{n — 1) 
2k 



ffl 



lim h 



m{n + 1) 
2k ' 



1 

m 



— 1 mh 

m 



m{n — 1) 
2k ' 



— 1 ffl(5_ 
m 



l/k 



5i ffl lim [£)fc-i/2(/im,„,fc)] 
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where 

is a centered random variable of variance n/m. In the regime when n = /3m, k = m? 
and m — )■ oo, the above measure converges weakly to a free difference of two identical 
Bernoulli random variables of parameter /3/2 and the result follows by the free central 
limit theorem. 

As for the formula (3) above, this is not exactly the random matrix result in [2], because 
the matrices studied there are of quite different nature of those studied here. We intend 
to come back to this phenomenon in some future work P|. 



5. Properties of the limiting measure 

In this section we study the atoms, support and density of the measures /im,n- These 
questions are of purely probabilistic nature, and our first result here is as follows. 

Proposition 5.1. The probability measure firn,n has at most one atom, at 0, of mass 
max{l — mn, 0}. 

Proof. The result follows from the general characterization of atoms of a free additive 
convolution, given by Bercovici and Voiculescu in a; is an atom for the free additive 
convolution of two measures z/i and 1^2 x = Xi + X2, where z/i({a;i}) + i'2{{x2}) > 1- 
And, in addition, if this is the case, then [ui ffl '^2]({a;}) = ^lii^i}) + ^2{{x2}) — 1- 

In our situation, since the free Poisson distribution of parameter c is known to have an 
atom at of mass max{l — c, 0}, we can apply the above-mentioned result when i>i is a 
free Poisson distribution of parameter m{n + l)/2 and 1^2 is the image of a free Poisson 
distribution of parameter m{n — l)/2 through the negation map, and we are done. □ 

Regarding now the support and density of firn,n, our main tool will be of course the 
degree 3 equation for the Cauchy transform appearing in Proposition 14. 2[ namely: 

- 1)(1 - G^) = mG{l + nC). 

In principle all the needed information can be obtained from these equations via the 
Stieltjes inversion formula, and various calculus methods for small degree polynomials. 

In what follows we will investigate a question which is particular interest for quantum 
information theory, as explained by Aubrun in P]: the problem is that of deciding, for 
various values of m > and n > 1, whether the support is contained in [0, 00) or not. 
The answer here is as follows (for a graphical representation, see Figure [3] below). 

Theorem 5.2. The measure ^m,n has positive support iff n < m/4 + 1/m and m > 2. 
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4 - 



3 - 



2 4 6 8 10 

Figure 3. Positivity of the support of fim,n- Measures with parameters 
belonging to the shaded region (m > 2, 1 < n < m/4 + 1/m) have positive 
support. 



Proof. We use the fact that the support of the a.c. part of the probabihty measure /im,n is 
a union of intervals defined by the points where the analyticity of G breaks. These points 
are the roots of the discriminant of the equation defining G (seen as an equation in G): 

A(0 = - 12m^^ + {n^m^ + Vlm^ - 20nm - 8)^^+ 
+ (-2n^m^ - 4m^ + 22nm^ - 20m)^+ 
+ n m — An m — 2nm + 12n m + m — 12nm + 4. 

This is a degree 4 equation in ^ that has either 2 or 4 real solutions (4 complex solutions 
is not a possibility since the support has to be non-empty). The discriminant of A(x) 
allows to decide between these cases: 

A2 (m,n) = -256m^(?2 - + 1) [m^n^ + ISm^n^ - 27m^ + 48mn - 64)^ . 

The sign of the discriminant of a quartic equation permits to decide whether the equation 
has two real and two complex roots or 4 roots of the same type (real or complex): the 
discriminant is negative iff the quartic has two real and two complex roots. In our case, 
the sign of A2 is the opposite of the sign of the bivariate polynomial 

P(m, n) = m^Tn? + 15m^n^ + 48mn — 27m^ — 64. 



Using straightforward calculus, it is easy to see that in the domain of interest here (m > 0, 
n > 1), P(m, n) < if and only if m < 4 and n < g{m) where g{m) is the only real root 
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1.5 r 




Figure 4. Behavior of the support of /im,n- Regions Al, A2 and B corre- 
spond to measures having 2 disjoint intervals for the a.c. support, whereas 
regions C and D have connected a.c. support. Measures with parameters 
in Al have a.c. support made of two intervals separated by the origin and 
measures in A2 have positive support. 



of P (seen as a polynomial in n): 



g[m) 



1 
2m 



/ X 1/3 / \ -1/3 

3.22/3 ^2 + m^ + mV4 + m2j + 6 • 2^/^ f 2 + + mV4 + m'^j -10 



The function g is larger than one for m G [0, 4] and it is decreasing in m in that range; its 
curve is plotted in Figure HI This settles the question of deciding if the support of /im,n 
has one or two disjoint intervals: measures corresponding to parameters in regions Al, 
A2 and B in Figure H] have support made of 2 disjoint intervals, whereas regions C and D 
correspond to measures having connected absolutely continuous support. 

Another indicator for positive support is the value at the origin ^ = of the discriminant 

A: 

A(0) = (m^ - 4mn + 4)(mn - 1)^. 

The two simple curves g{m) and h{m) = m/4 + 1/m partition the set (0, 00) x (1, 00) 
into five connected components, Al, A2, B, C, D, see Figure Si Points in the regions 
Al, A2 and D have A(0) > 0, whereas points in the regions B and C have A(0) < 0. 
A negative value at the origin indicates (since the leading coefficient of A(^) is strictly 
positive) that zero is in the interior of the support of /im,n- Hence, positively supported 
measures correspond to a subset of AlU A2U D. Moreover, since the average of firn,n is 
positive and measures in D have connected support, it follows that all the measures in D 
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Signs of A'(0), A"(0) 



A'(0) < 0, A"(0) < 
A'(0) < 0, A"(0) > 
A'(0) > 0, A"(0) < 
A'(0) > 0, A"(0) > 
Table 1. Nature of the roots of A'({ 



Nature of the roots of A' 



2 negative, 1 positive 

3 positive 

1 negative, 2 positive 
1 negative, 2 positive 
in terms of the signs of A'(0) and A"(0). 



are positively supported. The unique point of intersection of the graphs of g and h in the 
domain m>0, n>lisn = + ^ ^ 1.01372, m = {-23n + 3n^) ^ 2.35992. 

To conclude, one needs to understand how measures with parameters in Al and A2 
behave. We claim that measures in Al are supported on Ji U I2 where Ji C M_, I2 C M+ 
and measures in A2 are supported on Ji U I2 where Ji, /2 C M+. 

To show this, note that both regions Al and A2 are below g and h, hence the support has 
two connected components /i, I2 and it does not contain zero. Without losing generality, 
we assume that Ii is at the left of l2- Three cases need to be considered: zero is at the 
left of Ii, zero separates Ii and I2 or zero is at the right of J2. The last case cannot occur, 
since /i^.n cannot be supported entirely on the negative half line (it has positive mean). 
To distinguish between the first two cases, we analyze the derivatives of A, given by 

A'(^) = 16^^ - 36m^2 + {2n^m'^ + 24m^ - AOnm - 16) ^- 

- 2n^m^ - 4m^ + 22nm^ - 20m; 
A"(0 = 48^^ - 72m^ + 2n^m^ + 2^m^ - AOnm - 16. 

The sign of the roots of A' allows to distinguish between the two cases above: A has 4 
positive roots if and only if A' has 3 positive roots. The other case where A has two 
positive and two negative roots corresponds to A' having at least one (and at most two) 
negative roots. Since both A' and A" have negative sub leading coefficient, it follows that 
these polynomials must have at least one positive root. In Tabled! the nature of the roots 
of A' is given in terms of the sign of A'(0) and A"(0). 

Both equations A'(0) = and A"(0) = are of second degree in n, with solutions 



and respectively 



, , 11 ± V81 - 8m2 , , 

n = pi 2[m) = fn £ (0, 9/ v 8J 

2m 



10 ± 2V3V9^^ 

n = qio[m) = m G (0, 3 , 

m 



hence A'(0) > if and only if m G (0, 9/v^) and pi < m < p2 and A"(0) < if and only 
if m G (0, 3) and qi < m < q2- These functions are plotted in Figured 



BLOCK-TRANSPOSED WISHART MATRICES 



15 



n 




12 3 



Figure 5. Plots of the solutions for the equations A'(0) = and A"(0) = 
0, along with h{m). 

Straightforward computations show that pi{m) > 1 for m G (0, 5/3] U [2, 9/^8] and 
that h{m) > pi{m) for m G (0, 2). The equation gi(m) = 1 has a unique positive solution 
5/3 < 2(5 + \/5T)/13 < 2. From these facts, one can conclude that the only couples 
(m, n) for which A'(0) < and A"(0) > (or, equivalently, the equation A(^) = has 4 
positive real roots) and which are below the graph of h (i.e. satisfy n < m/A + 1/m) are 
those for which m > 2, which concludes the proof. □ 

6. Concluding remarks 

We have seen in this paper that Aubrun's recent and surprising result in [T] on the 
block-transposed Wishart matrices, which basically says that "when performing the block 
transposition the Marchenko-Pastur law becomes a shifted semicircle law" is best under- 
stood first by enlarging the set of parameters, and then by using free probability theory. 
Indeed, by replacing Aubrun's parameters {(P,a(P) by general parameters {dn,dm), and 
then by doing a free probabilistic study, our conclusion is that "when performing the 
block transposition, the free Poisson law becomes a free difference of free Poisson laws" . 

With this new point of view, several problems appear. First, since both in the usual 
Wishart and in the block-transposed Wishart cases we simply reach to certain "free linear 
combinations of free Poisson laws", one may wonder about a general result, covering 
both the Wishart and block-transposed Wishart cases. In addition, the random matrix 
result in [2], once again dealing with some "block- modified" Wishart matrices, and once 
again leading to certain free combinations of free Poisson laws, is waiting as well to be 
generalized. There are probably some connections as well with the work of Lenczewski 
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in [5]. We intend to come back to these algebraic questions, and to perform as well a 
systematic study of the resulting asymptotic measures, in a forthcoming paper [3]. 

As explained in the introduction, it is of interest for quantum information theory if 
the partial transpose of a typical matrix (or quantum state) is positive or not. For a 
fixed choice of parameters m and n, having an asymptotic eigenvalue counting measure 
which is supported on the positive half line is an indication that typical states are PPT. 
However, in order to decide if a quantum state is PPT or not, one needs to look at 
the smallest eigenvalue of the corresponding matrix. Aubrun has shown in [T] that the 
smallest eigenvalue converges indeed to the infimum of the support of the asymptotic 
measure; we conjecture that the same holds true in our setting. 
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